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Abstract 

We consider a partially asymmetric exclusion process (PASEP) on a finite num- 
ber of sites with open and directed boundary conditions. Its partition function was 
calculated by Blythe, Evans, Colaiori, and Essler. It is known to be a generating 
function of permutation tableaux by the combinatorial interpretation of Corteel and 
Williams. 

We prove bijectively two new combinatorial interpretations. The first one is 
in terms of weighted Motzkin paths called Laguerre histories and is obtained by 
refining a bijection of Foata and Zeilberger. Secondly we show that this partition 
function is the generating function of permutations with respect to right-to-left 
minima, right-to-left maxima, ascents, and 31-2 patterns, by refining a bijection of 
Frangon and Viennot. 

Then we give a new formula for the partition function which generalizes the 
one of Blythe & al. It is proved in two combinatorial ways. The first proof is 
an enumeration of lattice paths which are known to be a solution of the Matrix 
Ansatz of Derrida & al. The second proof relies on a previous enumeration of rook 
placements, which appear in the combinatorial interpretation of a related normal 
ordering problem. We also obtain a closed formula for the moments of Al-Salam- 
Chihara polynomials. 

1 Introduction 

1.1 The PASEP partition function 

The partially asymmetric simple exclusion process (also called PASEP) is a Markov chain 
describing the evolution of particles in sites arranged in a line, each site being either 
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empty or occupied by one particle. Particles may enter the leftmost site at a rate a > 0, 
go out the rightmost site at a rate (3 >0, hop left at a rate g > and hop right at a rate 
p > when possible. By rescaling time it is always possible to assume that the latter 
parameter is 1 without loss of generality. It is possible to define either a continuous-time 
model or a discrete-time model, but they are equivalent in the sense that their stationary 
distributions are the same. In this article we only study some combinatorial properties 
of the partition function. For precisions, background about the model, and much more, 
we refer to [5, 6, 11, 12, 16, 30]. We refer particularly to the long survey of Blythe and 
Evans [4] and all references therein to give evidence that this is a widely studied model. 
Indeed, it is quite rich and some important features are the various phase transitions, and 
spontaneous symmetry breaking for example, so that it is considered as a fundamental 
model of nonequilibrium statistical physics. 

A method to obtain the stationary distribution and the partition function of the 
model is the Matrix Ansatz of Derrida, Evans, Hakim and Pasquier [16]. We suppose that 
D and E are linear operators, {W\ is a vector, \V) is a linear form, such that: 



DE -qED = D + E, {W\aE = {W\, PD\V) = \V), {W\V) = 1, (1) 



then the non- normalized probability of each state can be obtained by taking the product 
{W\ti . . AnIV) where U is D if the zth site is occupied and E if it is empty. It follows that 
the normalization, or partition function, is given by (VF|(D + E)^\V). It is possible to 
introduce another variable y, which is not a parameter of the probabilistic model, but is 
a formal parameter such that the coefficient of in the partition function corresponds to 
the states with exactly k particles (physically it could be called a fugacity). The partition 
function is then: 



which we may take as a definition in the combinatorial point of view of this article (see 
Section 2 below for precisions). An interesting property is the symmetry: 



which can be seen on the physical point of view by exchanging the empty sites with 
occupied sites. It can also be obtained from the Matrix Ansatz by using the transposed 
matrices D* and E* and the transposed vectors {V\ and \W), which satisfies a similar 
Matrix Ansatz with a and (3 exchanged. 

In section 4, we will use an explicit solution of the Matrix Ansatz [5, 6, 16], and it will 
permit to make use of weighted lattice paths as in [6]. 

1.2 Combinatorial interpretations 

Corteel and Williams showed in [11, 12] that the stationary distribution of the PASEP 
(and consequently, the partition function) has a natural combinatorial interpretation in 
terms of permutation tableaux [32]. This can be done by showing that the two operators 



{W\{yD + Ef\V), 



(2) 



ZN{a,^,y,q) = y^ZN{f3,a, ^ g), 



(3) 
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D and E of the Matrix Ansatz describe a recursive construction of these objects. They 
have in particular: 

where PT^+i is the set of permutation tableaux of size + 1, a(T) is the number of Is 
in the first row, h{T) is the number of unrestricted rows, r{T) is the number of rows, and 
w{T) is the number of superfluous Is. See Definition 3.1.1 below, and [12, Theorem 3.1] 
for the original statement. Permutation tableaux are interesting because of their link 
with permutations, and it is possible to see Z^y as a generating function of permutations. 
Indeed thanks to the Steingn'msson- Williams bijection [32], it is also known that [12]: 

o-e6jv+i 

where we use the statistics in the following definition. 
Definition 1.2.1. Let a G ©„. Then: 

• u{a) the number of sj^ecza/ right-to-left minima, i.e. integers j G {l,...,n} such 
that a{i) = minj<i<„cr(i) and cr(j) < a{l), 

• v{a) is the number of special left-to-right maxima, i.e. integers j G {1, . . . ,n} such 
that a{j) = maxi<j<jCT(z) and a{j) > (t(1), 

• wex{a) is the number of weak exceedances of a, i.e. integers j G {1, . . . ,?7,} such 
that a{j)> j, 

• and cr(cr) is the number of crossings, i.e. pairs {i,j) G {1, . . . ,n}'^ such that either 
i<j< cr{i) < a{j) or a{i) < a{j) <i < j. 

It can already be seen that Stirling numbers and Eulerian numbers appear as special 
cases of Zj^. We will show that it is possible to follow the statistics in (5) through the 
weighted Motzkin paths called Laguerre histories (see [9, 33] and Definition 3.1.2 below), 
thanks to the bijection of Foata and Zeilberger [9, 19, 29]. But we need to study several 
subtle properties of the bijection to follow all four statistics. We obtain a combinatorial 
interpretation of in terms of Laguerre histories, see Theorem 3.2.4 below. Even more, 
we will show that the four statistics in Laguerre histories can be followed through the 
bijection of Frangon and Viennot [9, 20]. Consequently we will obtain in Theorem 3.3.3 
below a second new combinatorial interpretation: 

o-e6jv+i 

where we use the statistics in the next definition. This was already known in the case 
a = 1, see [9, 10]. 

Definition 1.2.2. Let a G 6„. Then: 
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• s(cr) is the number of right-to- left maxima of a and t{a) is the number of right-to- left 
minima of a, 

• asc((T) is the number of ascents, i.e. integers i such that either i = n or 1 < i < n — 1 

and a{i) < a{i + 1), 

• 31-2(cr) is the number of generalized patterns 31-2 in a, i.e. triples of integers 
(i, i + such that l<i<z + l< j<n and a{i + 1) < cr(j) < a{i). 



1.3 Exact formula for the partition function 

An exact formula for Zn was given by Blythe, Evans, Colaiori, Essler [5, Equation (57)] 
in the case where y = 1. It was obtained from the eigenvalues and eigenvectors of the 
operator D + E as defined in (16) and (17) below. This method gives an integral form for 
Zn, which can be simplified so as to obtain a finite sum rather than an integral. Moreover 
this expression for Zjy was used to obtain various properties of the large system size limit, 
such as phases diagrams and currents. Here we generalize this result since we also have 
the variable y, and the proofs are combinatorial. This is an important result since it is 
generally accepted that most interesting properties of a model can be derived from the 
partition function. 

Theorem 1.3.1. Let a = {I - q)- -1 and P = {1 - q)^ - 1. We have: 



'/3 

N 

i 



■jj-^--^^RN,n{y,q)Bnia,/3,y,q), (7) 



n=0 



where 



I N — n I 

L^— J N-n-2i 



3=0 



and 



Bn{a,l3,y,q) = ^ 



k=0 



n 



a\y^T-'. (9) 



In the case where y = 1, one sum can be simplified by the Vandermonde identity 

(^) (Jl •) = Cm) ' recover the expression given in [5, Equation (54)] by Blythe 
& al:^ 

RNAi,q)= ± (-i)^(U-l" J-UJ-V.))'z^"^)m,- (10) 

i=0 

In the case where a = /3 = 1, it was known [14, 23] that (1 — q)^^^Z]\f is equal to: 

N+l /N+l-k \ / ^ \ 

k=0 \ j=0 / \i=0 / 



4 



(see Remarks 4.3.3 and 5.0.6 for a comparison between this previous result and the new 
one in Theorem 1.3.1). And in the case where y = q = 1, from a recursive construction 
of permutation tableaux [10] or lattice paths combinatorics [6] it is known that: 



The first proof of (7) is a purely combinatorial enumeration of some weighted Motzkin 
paths defined below in (19), appearing from explicit representations of the operators D and 
E of the Matrix Ansatz. It partially relies on results of [14, 23] through Proposition 4.1.1 
below. In contrast, the second proof does not use a particular representation of the 
operators D and E, but only on the combinatorics of the normal ordering process. It also 
relies on previous results of [23] (through Proposition 5.0.4 below), but we will sketch a 
self-contained proof. 

This article is organized as follows. In Section 2 we recall known facts about the 
PASEP partition function Z^v, mainly to explain the Matrix Ansatz. In Section 3 we 
prove the two new combinatorial interpretations of Z^r, starting from (5) and using various 
properties of bijections of Foata and Zeilberger, Frangon and Viennot. Sections 4 and 5 
respectively contain the the two proofs of the exact formula for Zn in Equation (7). In 
Section 6 we show that the first proof of the exact formula for can be adapted to 
give a formula for the moments of Al-Salam-Chihara polynomials. Finally in Section 7 we 
review the numerous classical integer sequences which appear as specializations or limit 
cases of Z^- 
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2 Some known properties of the partition function 



As said in the introduction, the partition function Z^ can be derived by taking the product 
{W\{yD + E)^\V) provided the relations (1) are satisfied. It may seem non-obvious that 
{W\{yD + E)^\V) does not depend on a particular choice of the operators D and E, and 
the existence of such operators D and E is not clear. 

The fact that {W\{yD + E)^\V) is well-defined without making D and E explicit, 
in a consequence of the existence of normal forms. More precisely, via the commutation 
relation DE — qED = D + E we can derive polynomials c-^^ in y and q with non-negative 
integer coefficients such that we have the normal form: 




(12) 





(13) 
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(this is a finite sum). See [3] for otlier combinatorial interpretation of normal ordering 
problems. It turns out that the 4j ^ are uniquely defined if we require the previous equality 
to hold for any value of a, /3, y and q, considered as indeterminates. Then the partition 
function is: 

Z^ia,(3,y,q) = {W\iyD + Ef\V) = c'Hj' P'' ■ (14) 

Indeed, this expression is valid for any choice of \V), D and E since we only used 
the relations (1) to obtain it. In particular Zj^ is a polynomial in ?/, g, ^ and with 
non-negative coefficients. For convenience we also define: 

ZN{a,P,y,q) = ZN{^,^,y,q). (15) 

For example the first values are: 

Zq = 1, Zi = a + y(3, 2^ = 0?^ y{a + /3 + af3 + af3q) + y'^f3'^, 

Zg = y^P^ + {a(5'^q + a/?^ + a + a/3 + a^'^q^ + [5 + (3'^q + 2a(5q + 2(3^) 

+ {2a^ + a^g + a + /3a^g^ + + (3a'^q + a/3 + /3 + 2a/3q) y + a^. 

Even if it is not needed to compute the first values of Zn, it is useful to have explicit 
matrices D and E satisfying (1). The best we could hope is finite-dimensional matrices 
with non-negative entries, however this is known to be incompatible with the existence 
of phase transitions in the model (see section 2.3.3 in [4]). Let a = (1 — g)^ — 1 and 
/3 = (1 — g)-i — 1, a solution of the Matrix Ansatz (1) is given by the following matrices 
D = (A,i)i,jeN and E = {Eij)ij^n (see [16]): 



(1 - g) A,i = 1 + ^q\ 
(1 - q)Ei^i = 1 + aq\ 
all other coefficients being 0, and vectors: 

(1^1 = (1,0,0,...), 



(1 - g)A,i+i = 1 - ak\ (16) 
(l-g)A+i,, = l-g*+\ (17) 

|V) = (1,0,0,...)*, (18) 



{i.e. \V) is the transpose of {W\). Even if infinite-dimensional, they have the nice property 
of being tridiagonal and this lead to a combinatorial interpretation of Z^v in terms of lattice 
paths [6]. Indeed, we can see yD + E as & transfer matrix for walks in the non-negative 
integers, and obtain that (1 — q)'^ Zj^ is the sum of weights of Motzkin paths of length N 
with weights: 

• 1 — g'^"'"^ for a step starting at height /i, 

• (1 + 2/) + (a + yj3)q^ for a step — starting at height /i, (19) 

• y{l — ajSq^'"^) for a step \ starting at height h. 
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We recall that a Motzkin path is similar to a Dyck path except that there may be hori- 
zontal steps, see Figures 1, 3, 4, 5 further. These weighted Motzkin paths are our starting 
point to prove Theorem 1.3.1 in Section 4. 

We have sketched how the Motzkin paths appear as a combinatorial interpretation of 
Zn starting from the Matrix Ansatz. However it is also possible to obtain a direct link 
between the PASEP and the lattice paths, independently of the results of Derrida & al. 
This was done by Brak & al in [6], in the even more general context of the PASEP with 
five parameters. 

3 Combinatorial interpretations of Zn 

In this section we prove the two new combinatorial interpretation of Z^. Firstly we prove 
the one in terms of Laguerre histories (Theorem 3.2.4 below), by means of a bijection orig- 
inally given by Foata and Zeilberger. Secondly we prove the one in terms in permutations 
(Theorem 3.3.3 below). 

3.1 Permutation tableaux and Laguerre histories 

We recall here the definition of permutation tableaux and their statistics needed to state 
the previously known combinatorial interpretation (4). 

Definition 3.1.1 ([32]). Let A be a Young diagram (in English notation), possibly with 
empty rows but with no empty column. A complete filling of A with O's and I's is a 
permutation tableau if: 

• for any cell containing a 0, all cells above in the same column contain a 0, or all 
cells to the left in the same row contain a 0, 

• there is at least a 1 in each column. 

A cell containing a is restricted if there is a 1 above. A row is restricted if it contains 
a restricted 0, and unrestricted otherwise. A cell containing a 1 is essential if it is the 
topmost 1 of its column, otherwise it is superfluous. The size of such a permutation 
tableaux is the number of rows of A plus its number of columns. 

To prove our new combinatorial interpretations, we will give bijections linking the 
previously-known combinatorial interpretation (5), and the new ones. The main combi- 
natorial object we use are the Laguerre histories, defined below. 

Definition 3.1.2 ([33]). A Laguerre history of size n is a weighted Motzkin path of n 
steps such that: 

• the weight of a step starting at height h is yq'^ for some i E {0, . . . , h}, 

• the weight of a step — )■ starting at height h is either yq'- for some i E {0, . . . , h} or 

for some i E {0, . . . , h — 1}, 
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• the weight of a step \ starting at height h is for some i G {0, . . . , /i — 1}. 

The total weight of the Laguerre history is the product of the weights of its steps. We call 
a type 1 step, any step having weight yq^ where h is its starting height. We call a type 2 
step, any step having weight q^~^ where h is its starting height. 

As shown by P. Flajolet [18], the weighted Motzkin paths appear in various combi- 
natorial contexts in connexion with some continued fractions called J-fractions. We also 
recall an important fact from combinatorial theory of orthogonal polynomials. 

Proposition 3.1.3 (Flajolet [18], Viennot [33]). // an orthogonal sequence {PnjneN is 
defined by the three-term recurrence relation 

xPn{x) = Pn+l{x) + hnPn{x) + A„P„_i(x), (20) 

then the moment generating function has the J-fraction representation 

A? 

1 - felt - ^ 



equivalently the nth moment fi^ is the sum of weights of Motzkin paths of length n where 
the weight of a step (respectively — \j starting at height h is ah (respectively b^, Ch) 
provided A„ = an_ic„. 

Remark 3.1.4. The sum of weights of Laguerre histories of length n is the nth mo- 
ment of some g-Laguerre polynomials (see [25]), which are a special case of rescaled 
Al-Salam-Chihara polynomials. On the other hand is the A^th moment of shifted Al- 
Salam-Chihara polynomials (see Section 6). We will use the Laguerre histories to derive 
properties of Z^, however they are related with two different orthogonal sequences. 



3.2 The Foata-Zeilberger bijection 

Foata and Zeilberger gave a bijection between permutations and Laguerre histories in [19]. 
It has been extended by de Medicis and Viennot [29], and Corteel [9]. In particular, Corteel 
showed that through this bijection "^fz we can follow the number weak exceedances and 
crossings [9]. The bijection "^pz hnks permutations in 6„ and Laguerre histories of n 
steps. The ith step of "^pzicr) is: 

• a step ^ if i is a cycle valley, i.e. (J~^{i) > i < cr{i), 

• a step \ if z is a cycle peak, i.e. o-^^{i) < i > a'{i), 

• a step — )■ in all other cases. 



8 



And the weight of the ith step in "^pzicr) is y^q^ with: 



6 = 1 if i < a{i) and otherwise, 
j = if{k\k<i< a{k) < cr{i) } Hi < a{i), 
< a{k) <i <k} \i a{i) < i. 



#{ A: 



It follows that the total weight of "ifpzif^) is ywex(o-)gcr((T)_ ^j^g statistics wex and 

cr in a permutation a, it is practical to represent a by an arrow diagram. We draw n 
points in a line, and draw an arrow from the ith point to the cr(i)th point for any i. This 
arrow is above the axis if z < cr{i), below the axis otherwise. Then wex(cr) is the number 
of arrows above the axis, and cr((j) is the number of proper intersection between arrows 
plus the number of chained arrows going to the right. See Figure 1 for an example with 
cr = 672581493, so that wex(cr) = 5 and cr(cr) = 7. 












yq^ 










y<A 


















































Figure 1: The permutation a = 672581493 and its image '^fz{o'). 



Lemma 3.2.1. Let a G &n, and 1 < i < n. Then i is a left-to-right maximum of a if 
and only if the ith step o/\E'i?2((j) is a type 1 step (as in Definition 3.1.2). 

Proof. Let us call a (a, 'i)-sequence a strictly increasing maximal sequence of integers 
Ml, . . . , such that a{uj?) = lifc+i for any 1 < A; < j — 1, and also such that ui < i < Uj. 
By maximality of the sequence, Ui is a cycle valley and Uj is a cycle peak. The number 
of such sequences is the difference between the number of cycle valleys and cycle peaks 
among {1, . . . , i — 1}, so it is the starting height h of the zth step in '^fz{o'). 

Any left-to- right maximum is a weak exceedance, so i is a left-to-right maxima of a 
if and only if i < a{i) and there exists no j such that j < i < < cr{j). This is also 
equivalent to the fact that i < a{i), and there exists no two consecutive elements Uk, Uk+i 
of a (cr, i)-sequence such that Uk < i < < "Ufe+i. This is also equivalent to the fact 
that i < cr(z), and any (cr, z)-sequence contains two consecutive elements u^, Uk+i such 
that Uk < i < Uk+i < cr(i). 

By definition of the bijection "^pz it is equivalent to the fact that the zth step of 
^Fz(cr) has weight yq^, i.e. the ith step is a type 1 step. □ 

Lemma 3.2.2. Let o G ©„, and 1 < i < n. We suppose i ^ cr^i)- Then i is a right-to-left 
minima of a if and only if the ith step of "^pzicr) is a type 2 step. 
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Proof. We have to pay attention to the fact that a right-to-left minimum can be a fixed 
point and we only characterize the non-fixed points here. This excepted, the proof is 
similar to the one of the previous lemma. □ 



Before we can use the bijection \1/_fz we need a slight modification of the known 
combinatorial interpretation (5), given in the following lemma. 

Lemma 3.2.3. We have: 



where u'{a) is the number of right-to- left minima i of a satisfying a ^{N < i. 

Proof. This just means that in (5) we can replace the statistic u with u', and this can 
be done via a simple bijection. For any cr g &n+i, let a be the reverse complement of 
a~^, i.e. a{i) = j if and only if a{N + 2 — j) = N + 2 — i. It is routine to check that 
u{a) = u'{a), wex((T) = wex(5"), and v{a) = v{a). Moreover, one can check that the arrow 
diagram of a is obtained from the one of a by a vertical symmetry and arrow reversal, so 
that cr (cr) = cr((T). So (5) and the bijection a i— t- a prove (22). □ 

From Lemmas 3.2.1, 3.2.2, and 3.2.3 it possible to give a combinatorial interpretation 
of Zjv in terms of the Laguerre histories. We start from the statistics in & tv+i described 
in Definition 1.2.1, then from (22) and the properties of "^fz we obtain the following 
theorem. 

Theorem 3.2.4. The polynomial yZ]^ is the generating function of Laguerre histories of 
+ 1 steps, where: 

• the parameters y and q are given by the total weight of the path, 

• (5 counts the type 1 steps, except the first one, 

• a counts the type 2 steps which are to the right of any type 1 step. 

Proof. Let a G &n+i- The smallest left-to-right maximum of a is 1, and any other left-to- 
right maximum i is such that a{l) < <7{i). So 1 is the only left-to-right maximum which 
is not special. So by Lemma 3.2.1, v{a) is the number of type 1 steps in "^pzicr), minus 
1. 

Moreover, a^^{N is the largest left-to-right maximum of a. Let i be a right-to-left 
minimum of a such that a~^{N -\- 1) < i. We have i ^ o'li), otherwise a would stabilize 
the interval 1, . . . , A^-|- 1} and this would contradict a^^{N + 1) < i. So we can apply 
Lemma 3.2.2, and it comes that u'{a) is the number of type 2 steps in \E'^2(c"), which are 
to the right of any type 1 step. So (22) and the bijection "^fz prove the theorem. □ 




(22) 
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Before ending this subsection, we sketch how to recover a known result in the case 
q = from Theorem 3.2.4. This was given in Section 3.2 of [7] (see also Section 3.6 in 
[4]) and proved via generating functions. For any Dyck path D, let ret(D) be the number 
of returns to height 0, for example ret(/^\) = 1 and Tet{/^\/^\) = 2, and the empty 
path ■ satisfies ret(-) = 0. The result is the following. 

Proposition 3.2.5 (Brak, de Gier, Rittenberg). When y = 1 and q = 0, the partition 
function is Zn = '^{■^Y^^^^'^^^y'^^^^^^ where the sum is over pairs of Dyck paths (Di,D2) 
whose lengths sum to 2N. 

Proof. When q = Owe can remove any step with weight in the Laguerre histories. When 
y = 1, to distinguish the two kinds of horizontal steps we introduce another kind of paths. 
Let us call a bicolor Motzkin path, a Motzkin path with two kinds of horizontal steps •••> 
and — )■, and such that there is no > at height 0. From Theorem 3.2.4, if y = 1 and q = 
then /3Zi^ is the generating function of bicolor Motzkin paths M of length + 1, where: 

• there is a weight /3 on each step or — )■ starting at height 0, 

• there is a weight a on each step \ or • > starting at height 1 and being to the right 
of any step with a weight (3. 

There is a bijection between these bicolor Motzkin paths, and Dyck paths of length 2A^ + 2 
(see de Medicis and Viennot [29]). To obtain the Dyck path D, each step in the bicolor 
Motzkin path M is replaced with a sequence of two steps Similarly, each step — )■ 

is replaced with /^\, each step ■ > is replaced with each step \ is replaced with 

\\. When some step s G {/^, — \} in M has a weight (3 or a, and is transformed into 
steps (si, S2) G \}^ in D, we choose to put the weight /3 or a on Si. It appears 

that D is a Dyck path of length 2A^ + 2 such that: 

• there is a weight [3 on each step starting at height 0, 

• there is a weight a on each step \ starting at height 2 and being to the right of 
any step with weight 

Then D can be factorized into Di D2 \i where Di and D2 are Dyck paths whose 
lengths sum to 2N, and up to a factor /3 it can be seen that (3 (respectively a) counts 
the returns to height in Di (respectively D2). More precisely the /3s are on the steps 
starting at height but there are as many of them as the number of returns to height 0. 
See Figure 2 for a an example. □ 

3.3 The Frangon- Viennot bijection 

This bijection was given in [20]. We use here the definition of this bijection given in [9]. 
The map ^'i^v' is a bijection between permutations of size n and Laguerre histories of n 
steps. Let a G &n, j & ■ ■ ■ ,n} and k = <7{j). Then the kth step of \1/fv(o") is: 
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Figure 2: The bijection between M, D and (Di, D2)- 



• a step if A; is a valley, i.e. cr(j — 1) > o-(j) < cr(j + 1), 

• a step \ if A; is a peak, i.e. a{j — 1) < o"(j) > cr(j + 1), 

• a step — )■ if is a double ascent, i.e. a{j — 1) < cr(j) < a{j + 1), or a double descent, 
I.e. a(i-l) >o{j) >a(j + l). 

This is done with the convention that ain + 1) = n + 1, in particular n is always an 
ascent of a G Moreover the weight of the A;th step is y^cl where 5 = 1 if j is an 
ascent and otherwise, and i = 31-2(cr, j). This number 31-2(a", j) is the number of 
patterns 31-2 such that j correspond to the 2, i.e. integers i such that 1 < i + 1 < j and 
cr(z + 1) < o"(j) < cr{i). A consequence of the definition is that the total weight of \E'Fv(cr) 
•g ^asc(o-)^3i-2(o-)_ Figures 3 and 4 for examples. 
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1 2 3 4 5 6 7 



Figure 3: Example of the permutation 4371265 and its image by the Frangon-Viennot 
bijection. 

Lemma 3.3.1. Let a G 6„ and 1 < i < n. Then cr^^{i) is a right-to-left minimum of a 
if and only if the ith step of'^pvic) is a type 1 step. 

Proof. This could be done by combining the arguments of [20] and [9]. We sketch a proof 
introducing ideas that will be helpful for the next lemma. 

We suppose that j = (y^^{i) is a right-to-left minimum. So j is an ascent, and any v 
such that i > a{v) is such that v < j. The integer 31-2{a,j) is the number of maximal 
sequence of consecutive integers u,u -\- 1, . . . ,v such that a{u) > a{u + 1) > ■ ■ • > ), 
and a{u) > i > (t{v). Indeed, any of these sequences u, . . . ,v is such that v < j and 
so it is possible to find two consecutive elements k,k + 1 in the sequence such that 
a{k + 1) < a{j) < a{k), and these k,k -\- 1 only belong to one sequence. 
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We call a (cr, «)-sequence a maximal sequence of consecutive integers u,u + 1, . . . ,v 
such that a{u) > a{u + 1) > ■ ■ ■ > cr(f ), and a{u) > i > cr{v). By maximality, m is a peak 
and f is a valley. The number of such sequences is the difference between the number 
of peaks and number of valleys among the elements of image smaller than i, so it is the 
starting height h of the zth step in '^fv{o'). 

So with this definition, we can check that j = <y^^{i) is a right-to-left minimum of a 
if and only if j is an ascent and any (a, z)-sequence u,u + 1, ... ,v is such that v < j. So 
this is equivalent to the fact that the zth step of '^fv{<^) is a type 1 step. 



Lemma 3.3.2. Let a G (5„, and 1 < i < n. We suppose a ^{i) < n. Then a ^{i) is a 
right-to-left maximum of a if and only if 

• the ith step o/\E'fv'(o") it is a type 2 step, 

• any type 1 step is to the left of the ith step. 

Proof. We keep the definition of (a, z)-sequence as in the previous lemma. First we suppose 
that ^^^{i) is a right-to-left maximum strictly smaller than n, and we check that the two 
points are satisfied. If a~^{j) is a right-to-left minimum, then i > j, so the second point 
is satisfied. A right-to- left maximum is a descent, so the ith step is — t- or \ with weight 
q^. We have to show g = h — 1. Since cr~^(z) is a right-to-left maximum, there is no 
(cr, 2)-sequence u < ■ ■ ■ < v with <J~^{i) < u. So there is one (cr, i)-sequence u < ■ ■ ■ < v 
such that u < (T~^{i) < v, and the h — 1 other ones contains only integers strictly smaller 
than a~^{i). So the ith step of "i/Fvicr) has weight q'^'^. 

Reciprocally, we suppose that the two points above are satisfied. There are h — 1 (a, i)- 
sequence containing integers strictly smaller than a^^{i). Since cr~^(z) is a descent, the 
hth (cr, i)-sequence u < ■ ■ ■ < v is such that u < cr^^{i) < v. So there is no (a, i)-sequence 
u < ■ ■ ■ < V such that cr^^(i) < u. 

If we suppose that i is not a right-to-left maximum, there would exist k > i such that 
a^^{k) > cr^^(i). We take the minimal k satisfying this property. Then the images of 
cx^^{k) + 1, . . . , n are strictly greater than k, otherwise there would exist i > a^^{k) such 
that a{£) > i > a{£ + 1). But then a^^{k) would be a right-to-left minimum and this 
would contradict the second point that we assumed to be satisfied. □ 

In Theorem 3.2.4 we have seen that Z^- is a generating function of Laguerre histo- 
ries, and the bijection "^pv together with the two lemmas above give our second new 
combinatorial interpretation of Z^. 

Theorem 3.3.3. We have: 



□ 




,asc(o-)-l 31-2((t) 



(23) 



where we use the statistics in Definition 1.2.2 above. 
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For example, in Figure 4 we have a permutation a such that 

Indeed '^pv{H) has total weight y^q^ , has four type 1 steps and two type 2 steps to the 
right of the type 1 steps. 

9 
8 
7 
6 
5 
4 
3 
2 
1 

123456789 
Figure 4: The permutation a = 812563974 and its image by "^fv- 

Remark 3.3.4. We have mentioned in the introduction that the non-normalized prob- 
ability of a particular state of the PASEP is a product {W\ti . . .t]\[\V) . It is known 
[11] that in the combinatorial interpretation (4), this state of the PASEP corresponds to 
permutation tableaux of a given shape. It is also known [11] that in the combinatorial 
interpretation (5), this state of the PASEP corresponds to permutations with a given set 
of weak exceedances (namely, i -|- 1 is a weak exceedance if and only if ti = D, i.e. the ith 
site is occuppied). It is also possible to give such criterions for the new combinatorial in- 
terpretations of Theorems 3.2.4 and 3.3.3, by following the weak exceedances set through 
the bijections we have used. More precisely, in the first case the term {W\ti . . .tAr|F) 
is the generating function of Laguerre histories H such that ti = D if and only if the 
(A^ -|- 1 — z)th step in H is either a step — )■ with weight yq'' or a step \. In the second case, 
the term (VF|fi . . .tN\V) is the generating function of permutations a such that ti = D ii 
and only if a~^{N -|- 1 — i) is a double ascent or a peak. 

4 A first combinatorial derivation of Zn using lattice 
paths 

In this section, we give the first proof of Theorem 1.3.1. 

We consider the set *Piv of weighted Motzkin paths of length A^ such that: 

• the weight of a step starting at height h is q^ — g'^^ for some i E {0, . . . , h}, 

• the weight of a step — )■ starting at height h is either 1 + y or (a -|- yP)q^, 

• the weight of a step \ starting at height h is either y or —ya(3q^^^. 
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The sum of weights of elements in is (1 — (l)^ Z]^ because the weights sum to the 
ones in (19). We stress that on the combinatorial point of view, it will be important to 
distinguish (/i + 1) kinds of step starting at height /i, instead of one kind of step 
with weight 1 — q^^^. 

We will show that each element of *^Ar bijectively corresponds to a pair of weighted 
Motzkin paths. The first path (respectively, second path) belongs to a set whose generat- 
ing function is RN,n{y, q) (respectively, Bn{a, /3, y, q)) for some n G {0, . . . , N}. Following 
this scheme, our first combinatorial proof of (7) is a consequence of Propositions 4.1.1, 
4.2.1, and 4.3.1 below. 

4.1 The lattice paths for i?7v„(?/,g) 

Let DlN,n be the set of weighted Motzkin paths of length such that: 

• the weight of a step starting at height h is either 1 or — 

• the weight of a step — )■ starting at height h is either 1 + y or q^, 

• the weight of a step \ is y, 

• there are exactly n steps — )■ weighted by a power of q. 
In this subsection we prove the following: 

Proposition 4.1.1. The sum of weights of elements in ^N,n is R^^niViQ)- 

This can be obtained with the methods used in [14, 23], and the result is a consequence 
of the Lemmas 4.1.2, 4.1.3 and 4.1.4 below. 

Lemma 4.1.2. There is a weight-preserving bijection between ^N,n, o.nd the pairs {P,C) 
such that for some i G {0, . . . , ['^^y^J }, 

• P is a Motzkin prefix of length N and final height n + 2i, with a weight 1 + y on 
every step — and a weight y on every step \, 

• C is a Motzkin path of length n + 2i, such that 



the weight of a step starting at height h is 1 or —q^, 
the weight of a step — i- starting at height h is q^ , 
the weight of a step \ zs 1, 

there are exactly n steps — and no steps ^\ both with weights 1. 



(24) 



Proof. This is a direct adaptation of [14, Lemma 1]. 



□ 
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Lemma 4.1.3. The generating function of Motzkin prefixes of length N and final height 
n + 2i, with a weight 1 + y on every step — and a weight y on every step \, is 



N-n-2i 




Proof. This was given in [14, Proposition 4]. 



□ 



Lemma 4.1.4. The sum of weights of Motzkin paths of length n + 2i satisfying properties 



Some precisions are in order. In [14] and [23], we obtained the formula (11) which is 
the special case a = /3 = 1 in Z^r, and is the A^th moment of the g-Laguerre polynomials 
mentioned in Remark 3.1.4. Since Z^v is also very closely related with these polynomials 
(see Section 6) it is not surprising that some steps are in common between these previous 
results and the present ones. See also Remark 4.3.3 below. 

4.2 The lattice paths for Bn{a, f3,y,q) 

Let be the set of weighted Motzkin paths of length n such that: 

• the weight of a step starting at height h is either 1 or —q^^^, 

• the weight of a step — )■ starting at height h is {a + yP)q^, 

• the weight of a step \ starting at height h is —yaf3q'^~^. 
In this section we prove the following: 

Proposition 4.2.1. The sum of weights of elements in 23„ is Bn{a, (3,y,q). 

Proof. Let z/„ be the sum of weights of elements in 23^. It is homogeneous of degree n 
in a and /? since each step — )■ has degree 1 and each pair of steps and \ has degree 
2. By comparing the weights for paths in !B„, and the ones in (19), we see that Un is the 
term of (1 — q)"'Zn with highest degree in a and f3. Since a and (1 — g)^ (respectively, /? 
and (1 — q)^) only differ by a constant, it remains only to show that the term of Z„ with 

highest degree in a and /3 is Ylk=o [fc] Oi^{yP)^~'^- 

This follows from the combinatorial interpretation in Equation (4) in terms of permu- 
tation tableaux (see Definition 3.1.1). In the term of with highest degree in a and /3, 
the coefficient of a^j3^~^ is obtained by counting permutations tableaux of size n + 1, with 
n — k + 1 unrestricted rows. A; Is in the first row. Such permutation tableaux have n — k + 1 
rows, k columns, and contain no 0. They are in bijection with the Young diagrams that 
fit in a X (t2, — fc) box and give a factor [^] . □ 



(24) above is (-l)*g('^') [«+^] 



Proof. A bijective proof was given in [23, Lemmas 3, 4]. 



□ 
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We can give a second proof in relation with orthogonal polynomials. 

Proof. It is a consequence of properties of the Al-Salam-Carlitz orthogonal polynomials 
U^\x), defined by the recurrence [2, 27]: 

utlA^) = xUt\x) + (a + l)q'U^\x) + a{q' - ^^(x). (25) 

Indeed, from Proposition 3.1.3 the sum of weights of elements in 53„ is the nth moment 
of the orthogonal polynomial sequence {Pk{,x)}k>Q defined by 

Pfc+i(x) = xPk{x) + (a + yP)q''Pk + (g' - l)yah'''^Pk-i. (26) 

We have Pk{x) = {yPYu'^\x{yP)~^) where a = a{yf3)~^, and the nth moment of the 
sequence {uj^'^\x)}k>o is X]j=o [J] (^^^ §^ P]) ^^e article of D. Kim [26, Section 3] 
for a combinatorial proof). Then we can derive the moments of {Pk{x)}k>Oy and this gives 
a second proof of Proposition 4.2.1. □ 

Another possible proof would be to write the generating function 
continued fraction with the usual methods [18], use a limit case of identity (19.2.11a) in 
[15] to relate this generating function with a basic hypergeometric series and then expand 
the series. 

4.3 The decomposition of lattice paths 

Let ^H^„ be defined exactly as yiN,n, except that the possible weights of a step starting 
at height h are — g*"*"^ with i E {0, . . . , h}. The sum of weights of elements in is 
the same as with ^RN,n, because the possible weights of a step starting at height h sum 
to 1 — g^+^. Similarly let 03* be defined exactly as 23^, except that the possible weights 
of a step starting at height h are g* — g*"*"^ with z G {0, . . . , h}. 

Proposition 4.3.1. There exists a weight-preserving bijection $ between the disjoint 
union of ^Hjvn ^ over n G {0, . . . , A^}, and ^„ (we understand that the weight of a 
pair is the product of the weights of each element). 

To define the bijection, we start from a pair {Hi,H2) G x ^* for some n G 

{0, . . . , A} and build a path $(/7i, H2) G ^^n- Let z G {1, . . . , A}.' 

• If the ith step of Hi is a step — )■ weighted by a power of g, say the jth one among 
the n such steps, then: 

— the ith step $(i^i, H2) has the same direction as the jth step of H2, 

— its weight is the product of weights of the ith step of Hi and the jth step of 

H2. 

• Otherwise the ith step of $(i^i, H2) has the same direction and same weight as the 
ith step of Hi. 
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Figure 5: Example of paths Hi, H2 and their image H2) . 



See Figure 5 for an example, where the thick steps correspond to the ones in the first of 
the two cases considered above. It is immediate that the total weight of H2) is the 

product of the total weights of Hi and H2. 

The inverse bijection is not as simple. Let H G ^^n- The method consists in reading 
H step by step from right to left, and building two paths Hi and H2 step by step so 
that at the end we obtain a pair {Hi, H2) G x for some n e {0, . . . , A^}. At 

each intermediate stage, we have built two Motzkin suffixes, i.e. some paths similar to 
Motzkin paths except that the starting height may be non-zero. 

Let us fix some notation. Let if*--'^ be the Motzkin suffix obtained by taking the j 
last steps of H. Suppose that we have already read the j last steps of H, and built two 
Motzkin suffixes Hi'' and -^2''^ We describe how to iteratively obtain Hi^^^ and h'^^^'' ■ 
Note that H^'' and H^"^ are empty paths. Let h, h', and h" be the respective starting 
heights of H^^\ h[^^ and H^^\ 

This iterative construction will have the following properties, as will be immediate 
from the definition below: 

• ifp^ has length j, and the length of H2^ is the number of steps — t- weighted by a 
power of q in H^\ 

• We have h = h' + h". 

• The map $ as we described it can also be defined in the same way for Motzkin 
suffixes, and is such that H^^^ = ^{h[^\h[^^). 

To obtain H^^^"^ and H2^^\ we read the {j + l)th step in H starting from the right, and 
add steps to the left of H^^ and H2^ according to the following rules: 
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step read in H 


step added to i/^^ 


step added to 




— )• q^ 


\ -yoipq 


\ y 


\ y 




l + y 


l + y 




{a -\- yp)q 


h' 

— >■ q 


— > l^a + yp}q 


/• q' - with i < h! 






/ q' - q'^^ with i > h' 


h' 


gi—h' gi+l—h' 



We can also iteratively check the following points. 

• With this construction i/j''^^^ and H2^^^ are indeed Motzkin suffixes. This is be- 
cause we add a step to H^'' only in the case where i < h' , hence h' > 0. And 
we add a step to only in the case where i > h', hence h" > (since 
h = h' + h" > i). 

• The paths ifp^^^ and H2^^^ are respectively suffixes of an element in ^N,n and 
for some n G {0, . . . , A^}, i.e. the weights are valid. 

• The set of rules we have given is the only possible one such that for any j we have 
H(J) = ^h['\h^2'^). 

It follows that {Hl^\Hi'^^) G x 55^ for some n G {0, . . . , A^}, these paths are such 

that ^{h[^\h^2^^) = H, and it is the only pair of paths satisfying these properties. There 
are details to check, but we have a full description of $ and of the inverse map See 
Figure 6 for an example of the Motzkin suffixes we consider. 




Figure 6: Example of Motzkin suffixes used to define $ 



Before ending this subsection we can mention another argument to show that and 
the disjoint union of ^jvn ^ ®n have the same cardinal. Thus we could just focus on the 
surjectivity of the map $ and avoid making the inverse map explicit. The argument uses 
the notion of histories [33] and their link with classical combinatorial objects, as we have 
seen in the previous section with Laguerre histories. As an unweighted set, is a set of 
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colored Motzkin paths, with two possible colors on the steps — )■ or \, and h + 1 possible 
colors for a step starting at height h. So ^Pa? is in bijection with colored involutions / 
on the set {1, . . . , A^}, such that there are two possible colors on each fixed point or each 
arch (orbit of size 2). So they are also in bijection with pairs (/i,/2) such that for some 
ne{0,...,iV}: 

• Ii is an involution on {1, . . . , A^} with two possible colors on the fixed points (say, 
blue and red), and having exactly n red fixed points, 

• I2 is an involution on {1, . . . , n}. 

Using histories again, we see that the number of such pairs {h^h) is the cardinal of 

Remark 4.3.2. Note that considering ^P^r as an unweighted set is not equivalent to 
setting the various parameters to 1. For example the two possible colors for the horizontal 
steps correspond to the possible weights 1 + ?/ or (a + y[3)q\ This bijection using colored 
involution is not weight-preserving but it might be possible to have a weight-preserving 
version of it for some adequate statistics on the colored involutions. 

Remark 4.3.3. The decomposition $ is the key step in our first proof of Theorem 1.3.1. 
This makes the proof quite different from the one in the case a = (3 = 1 [14], even though 
we have used results from [14] to prove an intermediate step (namely Proposition 4.1.1). 
Actually it might be possible to have a direct adaptation of the case a = P = 1 [14] to 
prove Theorem 1.3.1, but it should give rise to many computational steps. In contrast 
our decomposition $ explains the formula for Zn as a sum of products. 



5 A second derivation of using the Matrix Ansatz 

In this section we build on our previous work [23] to give a second proof of (7). In this 
reference we define the operators 

D = + -/ and E = '^^E + -/, (27) 

q q q q 

where / is the identity. Some immediate consequences are 

DE-qED = ^—-^, {W\E = --{W\, and D\V) = --\V), (28) 
q'^ q q 

where a and /3 are defined as in the previous sections. While the normal ordering problem 
for D and E leads to permutation tableaux, for D and E it leads to rook placements as 
was shown for example in [35] . The combinatorics of rook placements lead to the following 
proposition. 
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Proposition 5.0.4. We have: 
(WHgyb + ,,E)''\V} 



where 



M=0 



i+j<k 
i+j=k mod 2 



'k-2l + u 
u 



i + j 
i 



(29) 



u 



k 

i-u-1 



(30) 



Proof. This is a consequence of results in [23] (see Section 2, Corollary 1, Proposition 12). 
We also give here a self-contained recursive proof. We write the normal form of {yqD + 
qE)'' as: 

{yqb + qEf = ^ d!^^] {qE)\qyby . (31) 
From the commutation relation in (28) we obtain: 



iqybyiqE) = q^{qE){qyby + y{l - q^){qyb) 



(32) 



If we multiply (31) by yqb + qE to the right, using (32) we can get a recurrence relation 



for the coefficients (i^-j , which reads: 



(fc+i) _ ,(fc) 



The initial case is that d-^ is 1 if (i, j) = (0, 0) and otherwise. It can be directly checked 
that the recurrence is solved by: 



(fc) 



i + j 
i 



M k-i-j 



,k 



(34) 



where we understand that M k~i~j jr. is when k — i — j is not even. More precisely, if we 



let e) 



(fc) 



k-i-3 



then we have: 



and also 



y{l~q^+')e'^^_,, = y{l-q 



1 + ] 
i 



M k-i-j + 1 



i + j 
i 



M k-i-j-l 



k- 



So to prove df'j = ef^j it remains only to check that 



(35) 



(36) 



(37) 



See for example [23, Proposition 12] (actually this recurrence already appeared more than 
fifty years ago in the work of Touchard, see loc. cit. for precisions). □ 
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Now we can give our second proof of Theorem 1.3.1. 
Proof. From (2) and (27) we have that (1 — q)^ Zj^ is equal to 

^ /AT\ 

[l + y)^-\-lf{W\{qyb + qEf\V). 



{Wm + y)I - qyb - qEf\V) = g 



So, from Proposition 5.0.4 we have: 



N 



fc=0 i+j<k 

i+j=k mod 2 



i + j 
i 



:i + y)''-''M,_^ 



(the (—1)^ cancels with a (—I)*"'"-'). Setting ri = i + we have: 



N 



ra=0 n<k<N ^ ^ 



n<k<N 
k=n mod 2 



So it remains only to show that the latter sum is RN,n{y, q)- If we change the indices so 
that k becomes n + 2/c, this sum is: 



fc=0 

N-n I 



N-n-2kyk ^(-1' 
i=0 



^qir) 



n + i 
i 



- C)) 



j=0 



n + i 
i 



We can simplify the latter sum by Lemma 5.0.5 below and obtain RN^n{y,q)- This com- 
pletes the proof. □ 

Lemma 5.0.5. For any N,n,i>0 we have: 



E 

k=i 



y 



k—i 



N 
n + 2k 



N-n~2k ( (n+2k 



N-n-2i 



(38) 



E/ ((j ) in+2i+j) ij-l) (n+2j+j+l)) 

j=0 



Proof. As said in Lemma 4.1.3, the right-hand side of (38) is the number of Motzkin 
prefixes of length A^, final height n + 2i, and a weight 1 + y on each step — )■ and y on each 
step \. Similarly, y^~^{{"'jj^^) — (^"_^?'\)) is the number of Dyck prefixes of length n + 2k 
and final height n + 2i, with a weight y on each step \. From these two combinatorial 
interpretations it is straightforward to obtain a bijective proof of (38). Each Motzkin 
prefix is built from a shorter Dyck prefix with the same final height, by choosing where 
are the N — n — 2k steps □ 
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Remark 5.0.6. All the ideas in this proof were present in [23] where we obtained the 
case a = j3 = 1. The particular case was actually more difficult to prove because several 
g-binomial and binomial simplifications were needed. In particular, it is natural to ask if 
the formula in (11) for ZAr|a=^=i can be recovered from the general expression in Theo- 
rem 1.3.1, and the (affirmative) answer is essentially given in [23] (see also Subsection 6.2 
below for a very similar simplification). 



6 Moments of Al-Salam-Chihara polynomials 

The link between the PASEP and Al-Salam-Chihara orthogonal polynomials Qn{x] a,b \ q) 
was described in [30]. These polynomials, denoted by Qn{x) when we don't need to specify 
the other parameters, are defined by the recurrence [27]: 

2xQ„(x) = Qn+i{x) + (a + b)q^Q^{x) + (1 - g")(l - abq^-')Q^^,{x) (39) 

together with Q-i{x) = and Qo{x) = 1. They are the most general orthogonal sequence 
that is a convolution of two orthogonal sequences [1]. They are obtained from Askey- 
Wilson polynomials Pn{x; a, b,c,d \ q) by setting c = = [27]. 



6.1 Closed formulas for the moments 

Let Qn{,x) = Qn(f — 1; tt, /? I q), where a = (1 — g)^ — 1 and /3 = (1 — g)-| — 1 as before. 

From now on we suppose that a = a and b = (5 (note that a and b are generic if a and j3 
are). The recurrence for these shifted polynomials is: 

xQn{x) = g„+i(x) + (2 + ag" + /3g")Q„(x) + (1 - g")(l - ah^-')Qn^r{x). (40) 

From Proposition 3.1.3, the A^th moment of the orthogonal sequence {Qn{x)}n>o is the 
specialization of (1 — q)^ Zn at y = 1. The A^th moment fiN of the Al-Salam-Chihara 
polynomials can now be obtained via the relation: 



N / \ 
I — n V / 



fc=0 

Actually the methods of Section 4 also give a direct proof of the following. 
Theorem 6.1.1. The Nth moment of the Al-Salam-Chihara polynomials is: 



0<n<Af \ 1=0 ^ ^ 



-J 



n=Nmod 2 ^ ' (41) 



X 
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Proof. The general idea is to adapt the proof of Theorem 1.3.1 in Section 4. Let C ^Pa? 
be the subset of paths which contain no step — )■ with weight 1 + y. By Proposition 3.1.3, 
the sum of weights of elements in ^'^y specialized at y = 1, gives the A^th moment of the 
sequence {Qn(f )}n>o- This can be seen by comparing the weights in the Motzkin paths 
and the recurrence (39). But the A'^th moment of this sequence is also 2^ fi^. 

From the definition of the bijection $ in Section 4, we see that ^{Hi, H2) has no step 
— > with weight 1 + y if and only if Hi has the same property. So from Proposition 4.3.1 the 
bijection gives a weight-preserving bijection between and the disjoint union of 
9^'^^ X 03* over n G {0, . . . , A^}, where C is the subset of paths which contain 
no horizontal step with weight 1 + y. Note that 5^7Vn is empty when n and don't have 
the same parity, because now n has to be the number of steps — )■ in a Motzkin path of 
length A^. In particular we can restrict the sum over n to the case n = N mod 2. 

At this point it remains only to adapt the proof of Proposition 4.1.1 to compute the 
sum of weights of elements in ^H'TVn) ^i^d obtain the sum over j in (41). As in the previous 
case we use Lemma 4.1.2 and Lemma 4.1.4. But in this case instead of Motzkin prefixes 
we get Dyck prefixes, so to conclude we need to know that ((^_.^/2_ ) ~ ((Ar_„w2_ _i) is 



the number of Dyck prefixes of length A^ and final height n + 2i. The rest of the proof is 



We have to mention that there are analytical methods to obtain the moments /xat of 
these polynomials. A nice formula for the Askey- Wilson moments was given by Stanton 
[31], as a consequence of joint results with Ismail [22, equation (1.16)]. As a particular 
case they have the Al-Salam-Chihara moments: 



where we use the g-Pochhammer symbol. The latter formula has no apparent symmetry 
in a and b and has denominators, but Stanton [31] gave evidence that (42) can be simpli- 
fied down to (41) using binomial, g-binomial, and g-Vandermonde summation theorems. 
Moreover (42) is equivalent to a formula for rescaled polynomials given in [25] (Section 4, 
Theorem 1 and equation (29)). 

6.2 Some particular cases of Al-Salam-Chihara moments 

When a = f) = in (41) we immediately recover the known result for the continuous q- 
Hermite moments. This is if A^ is odd, and the Touchard-Riordan formula if A^ is even. 
Other interesting cases are the g-secant numbers E2n{q) and g-tangent numbers E2n+i{q), 




similar. 



□ 




(42) 
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defined in [21] by continued fraction expansions of the ordinary generating functions: 

1 ^ , , 1 



n>0 



1 



and ^E2n+i{qY 



n>0 



(43) 



[2] t 



1 - 



The exponential generating function of the numbers i?n(l) is the function tan(a;) + sec(a;). 
We have the combinatorial interpretation [21, 24]: 

En{q) = Yl ^''"'^'^^ (44) 

where 2t„ C ©„ is the set of alternating permutations, i.e. permutations a such that 
cr(l) > cr(2) < (j(3) > . . . . The continued fractions show that these numbers are particular 
cases of Al-Salam-Chihara moments: 

E2n{q) = (i3^)^>2n|a=-6=i^, and £^2n+l(g) = ij^f''fJ'2n\a=-b=ig (45) 

(where i"^ = —1). From (41) and a g-binomial identity it is possible to obtain the closed 
formulas for i?2n('?) and E2n+i{q) that were given in [24], in a similar manner that (7) can 
be simplified into (11) when a = j3 = 1. Indeed, from (41) we can rewrite: 



a 



2m- 2 j 



(46) 



m=0 j,fc>0 

This latter sum over j and k is also 



E(-i)''''"''b:71P"1(^)' 



a 



2m- 2 j 



EM 



-l)'9''")p"r],L'],(;)'"'«'"-'M47) 

j,fc>o ^>i>o 

The sum over j can be simplified in the case a = —b = i^., or a = —b = iq, using the 
g-binomial identities already used in [23] (see Lemma 2): 



2m — j 




'£ 


£ 


q 


J. 



£(2m-£) 



(4J 



and 



Omitting details, this gives a new proof of the Touchard-Riordan-like formulas [24]: 



^ n 2m 



^(2m-£)+m 



m=0 



and 



£=0 
2m+l 



m=0 £=0 



(50) 



(51) 
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7 Some classical integer sequences related to 

It should be clear from the interpretation given in (6) that the polynomial Z^v contains 
quite a lot a of combinatorial information. When a = (3 = 1, the coefficients of in Zn 
are the g-Eulerian numbers introduced by Williams [36]: 

N 

ZN\a=i3=i = ^^y^Ek+i^N+iiq), (52) 

k=0 

where Ek,n{<l) is defined in [36, Section 6]. It was proved by Williams, that Ek,n{<l) is 
equal to the Eulerian number An^k when g = 1, to the binomial coefficient when 
q = —1, and to the Narayana number Nn± = ^{^){k-i) when q = 0. With the other 
parameters a and /3, there are other interesting results. 

7.1 Stirling numbers 

Carlitz g-analog of the Stirling numbers of the second kind, denoted by S2[n, k], are defined 
when 1 < k < nhj the recurrence [8]: 

S2[n, k] = S2[n — 1, k — 1] + [k]gS2[n — 1, A;], S'2[n, /c] = 1 if A; = 1 or A; = n. (53) 

Proposition 7.1.1. If a = 1, the coefficient of P^y^ in Z]^ is S'2[A^ + 1, A; + 1] . 

Proof. This follows from the interpretation (4) in terms of permutation tableaux (see 
Definition 3.1.1). Indeed, the coefficient of P^y^ in Zn counts permutation tableaux of 
size + 1, with A; + 1 rows, and A; + 1 unrestricted rows. In a permutation tableau with 
no restricted row, each column contains a sequence of O's followed by a sequence of I's. 
Such permutation tableaux follow the recurrence (53) where n is the size and k is the 
number of rows. Indeed, if the bottom row is of size we can remove it and this gives 
the term 5*2 [n — 1 , A; — 1] . Otherwise the first column is of size A;, this gives the term 
[A;]gS'2[n — 1, A;] because the factor [A;]^ accounts for the possibilities of the first column, 
the factor 5*2 [ti — 1, A;] accounts for what remains after removing the first column. □ 

The proof only relies on simple facts about permutation tableaux, and with no doubts 
it was previously noticed that 5*2 [n, A;] appears when we count permutation tableaux with- 
out restricted rows. Actually permutation tableaux with no restricted rows are equivalent 
to the 0-1 tableaux introduced by Leroux [28] as a combinatorial interpretation of 5*2 [n. A;]. 

From (7), it is possible to obtain a formula for 5*2 [^i-, A;]. First, observe that the coeffi- 
cient of y^ in Zn has degree k in /3. Hence, from the previous proposition: 

N 

y^a''S2[N + l,k + l] = hm Z^(l, 2 2/, g). (54) 

fe=0 
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We have i?7v,n(0, q) = (^) . When a = 1 and /3 = |, we have a = —q and y/3 = {l — q)a+y. 



So from (7) and (54) it is straightforward to obtain: 

N-k 



S2[N + l,k + l] 



1 



N-k 



3=0 





'k + f 


+ J J 


J 



(55) 



Note that this differs from the expression originally given by Carlitz [8]: 



S2[N,k] 



N-k 



N-k 
j=0 





'k + ] 




3 



(56) 



-I q 



but it is elementary to check that (55) and (56) are equivalent, using the two-term recur- 
rence relations for binomial and g-binomial coefficients. 

When y = a = 1, the coefficient of (3^ in is a g-analog of the Stirling number 
of the first kind 5*1 (A^ + -|- 1). It is such that q counts the number of patterns 
31-2 in permutations of size + 1 and with k + 1 right-to-left minima. Knowing the 
symmetry (3), we could also say that it is such that q counts the number of patterns 31-2 
in permutations of size -|- 1 and with A; + 1 right-to-left maxima. The combinatorial 
way to see the symmetry (3) is the transposition of permutation tableaux [13], so at the 
moment it is quite indirect to see that the two interpretations agree since we need all 
the bijections from Section 3. We have no knowledge of previous work concerning these 
g-Stirling numbers of the first kind. 



7.2 Fine numbers 

The sequence of Fine numbers shares many properties with the Catalan numbers, we 
refer to [17] for history and facts about them. We will show that a natural symmetric 
refinement of them appears as a specialization of Z^. 

A peak of a Dyck path is a factor /'\, we denote by pk(P) the number of peaks of a 
path P. A Fine path is a Dyck path D such that there is no factorization D = Di /^'\ D2 
where Di and D2 are Dyck paths. The Fine number F„ is the number of Fine path of 
length 2n, and more generally the polynomial Fn{y) is '^y^^^^^ where the sum is over 
Fine paths P of length 2n. These polynomials were considered in [17] via their generating 
function. 

An interesting property is that Fniy) is self- reciprocal, i.e. F^iy) = y"F„(i) (a simple 
proof of this will appear below). This is reminiscent of the Dyck paths: the number 
of Dyck paths of length 2n with k peaks is the Narayana number N^^^ and we have 
Nk,n = Nn+i-k,n- The first values are: 

Fi(y) = 0, F2iy)=y, F^{y) = y^ + y, F^{y) = y^ + Ay^ + y , 
F,{y) = y^ + ^y' + ^y'' + y. F^{y) = y'> + l?>y^ + 2<dy^ + l?>y^ + y. 

Proposition 7.2.1. When - = —y, g = 0, and (3 = 1, we have Zfq = F^iy). 
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Proof. In this case we have 13 = 0, a = —1 — y. From the weights in the general case (19), 
we see that now Zj\f is the sum of weights of Motzkin paths such that: 

• the weight of a step is 1, the weight of a step \ is y, 

• the weight of a step — )■ is 1 + but there is no such step at height 0. 

Let H{t,y) = 'EN>o^Nt^- I* is such that H{t,y) = I + yt^G{t,y)H{t,y), where G{t,y) 
counts the paths with the same weights but possibly with steps — )• at height 0. Let 
L{t,y) = yP^(^^ where i{P) is half the number of steps of P, and the sum is 

over all Dyck paths P. Some standards arguments show that G{t, y) is linked with 
Narayana numbers in such a way that L{t,y) = 1 + ytG{t,y). So we have H{t,y) = 
1 + f y) — l)H{t, y), which is precisely the functional equation given in [17, Section 7] 
for the generating function 'Y^Fn{y)t^- This completes the proof. □ 

When we substitute y with - in the Motzkin paths considered in the proof, we see 
that the weight of a step — )■ is divided by y and the weight of a step \ is divided by y"^, 
so the total weight is divided by where n is the length of the path. This proves the 
symmetry of the coefficients of Fn{y). 

Note that the symmetry of obtained in this section is not a particular case of pre- 
viously known symmetry (3). It may be a special case of another more general symmetry. 

8 Concluding remarks 

We have used two kinds of weighted Motzkin paths to study Z^. The first kind are the 
elements of ^^n, i-^- the paths coming from the matrices D and E defined in (16) and 
(17). They have the property that the weight of a step only depends on its direction and 
its height, so that there is a J- fraction expansion for the generating function 'Y1,n>o ^Nt^ 
with the four parameters a, (3, y and q. The second kind of weighted Motzkin paths 
are the Laguerre histories, and their nice property is that they are linked bijectively with 
permutations. One might ask if there is a set of weighted Motzkin paths having both 
properties, but its existence is doubtful. Still it could be nice to have a direct simple 
proof that these two kinds of paths give the same quantity Z^. 

Our two new combinatorial interpretations in Theorems 3.2.4 and 3.3.3 complete the 
known combinatorial interpretations (4) and (5), and this makes at least four of them. 
Although all is proved bijectively, there is not a direct bijection for any pair of combi- 
natorial interpretations. In particular it would be nice to have a more direct bijection 
between permutation tableaux and permutations to link the right-hand sides of (4) and 
(6), instead of composing four bijections (Steingn'msson- Williams, reverse complement of 
inverse, Foata-Zeilberger and Frangon-Viennot). Permutation tableaux are mainly inter- 
esting because of their link with permutations, so in this regard it is desirable to have a 
direct bijection preserving the four parameters considered here. 

We have given evidence that the lattice paths are good combinatorial objects to study 
the PASEP with three parameters. However, our combinatorial interpretation of Z^ with 
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the Laguerre histories rehes on the previous one with permutation tableaux. To complete 
the lattice paths approach, it might be interesting to have a direct derivation of stationary 
probabilities in terms of Laguerre histories. For example in [11], Corteel and Williams 
define a Markov chain on permutation tableaux which projects to the PASEP, similarly 
we could hope that there is an explicit simple description of such a Markov chain on 
Laguerre histories. 

The three-parameter PASEP is now quite well understood since we have exact expres- 
sions for many interesting quantities. In a more general model, we allow particles to enter 
the rightmost site, and exit the leftmost site, so that there are five parameters. In this 
case the partition function is linked with the Askey- Wilson moments, in a similar manner 
that the three-parameter partition function is linked with Al-Salam-Chihara moments 
[34]. Recently, Corteel and Williams [13] showed that there exist some staircase tableaux 
generalizing permutation tableaux, arising from this general model with five parameters. 
It is not clear whether a closed formula for the five-parameter partition function exists, 
and in the case it exists it might be unreasonably long. But knowing the results about 
the three-parameter partition function, we expect the five-parameter partition function 
to be quite full of combinatorial meaning. 
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